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ABSTRACT 

We postulate that supermassive black-holes grow in the centers of galaxies until they unbind the 
galactic gas that feeds them. We show that the corresponding self-regulation condition yields a correlation 
between black-hole mass (A/bh) and galaxy velocity dispersion (cr) as inferred in the local universe, and 
recovers the observed optical and X-ray luminosity functions of quasars at redshifts up to z ~ 6 based 
on the hierarchical evolution of galaxy halos in a ACDM cosmology. With only one free parameter and 
a simple algorithm, our model yields the observed evolution in the number density of optically bright or 
X-ray faint quasars between 2 < z < 6 across 3 orders of magnitude in bolometric luminosity and 3 orders 
of magnitude in comoving density per logarithm of luminosity. The self-regulation condition identifies 
the dynamical time of galactic disks during the epoch of peak quasar activity {z ~ 2.5) as the origin of 
the inferred characteristic quasar lifetime of ~ 10^ years. Since the lifetime becomes comparable to the 
Salpeter e-folding time at this epoch, the model also implies that the Mbh — c relation is a product of 
feedback regulated accretion during the peak of quasar activity. The mass-density in black-holes accreted 
by that time is consistent with the local black-hole mass density phh ^ (2-3 j^^ ^) x lO^M0Mpc^^, which we 
have computed by combining the Mbh-cr relation with the measured velocity dispersion function of SDSS 
galaxies (Sheth et al. 2003). Comparison of the local black-hole mass- function with that inferred from 
combining the feedback-relation with the halo mass-function suggests that most massive (> IO^Mq) 
black-holes may have already been in place by z ~ 6. Applying a similar self- regulation principle to 
supernova-driven winds from starbursts, we find that the ratio between the black hole mass and the 
stellar mass of galactic spheroids increases with redshift as (1 -I- z)'^/^ although the Mbh-c relation is 
redshift-independent . 

Subject headings: black-hole physics - quasars: general 



1. INTRODUCTION 

It has long been recognized that the formation processes 
of supermassive black-holes (BHs) and stars may be self- 
regulating (Dekel & Silk 1986; Silk & Rees 1998), since 
the amount of energy released by quasars and stars can ex- 
ceed the binding energy of the gas clouds in which they are 
born. Indeed, the stellar and central BH masses of galaxies 
show striking correlations with the depth of the gravita- 
tional potential well of their hosts. These phenomenolog- 
ical correlations take the form of the mass-velocity dis- 
persion correlation for BHs (Merritt & Ferrarese 2001; 
Tremaine et al. 2002) and the TuUy-Fisher (1977) and 
Faber-Jackson (1976) relations for stars in the present-day 
universe. 

In this paper we explore the self-limiting growth of the 
BH and stellar masses of galaxies as a function of cos- 
mic time, and investigate the consequences for the quasar 
luminosity function. Recent work suggests that the ki- 
netic output of quasars may be comparable to their radia- 
tive output (see §2.1.1 in Furlanetto & Loeb 2001; Begel- 
man 2003). Our simple approach postulates that growth 
stops when the energy released by these sources exceeds 
the binding energy of the gas within their host galaxies 
after one dynamical time of the gas reservoir. Using the 
redshift-dependent abundance of galactic halos as a func- 
tion of potential depth in a ACDM cosmology, we demon- 
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strate that this simple approach describes the luminosity 
function of quasars in optical and X-ray bands. We adopt 
the WMAP set of cosmological parameters (Bennett et al. 
2003), namely mass density parameters of 17^, = 0.27 in 
matter, J7b = 0.044 in baryons, JIa = 0.73 in a cosmologi- 
cal constant, a Hubble constant of Hq — 71 kms~^ Mpc^^, 
an rms amplitude of erg = 0.84 for mass density fluctua- 
tions in a sphere of radius 8/i~^Mpc, and a primordial 
power-spectrum with a power- law index n = 1. 

We begin the paper with a discussion of the regulation 
of BH growth by feedback during luminous quasar phases 
(§ 2) and then discuss the implications for modeling of the 
optical (§ 3.1) and X-ray (§ 3.2) luminosity functions of 
quasars, including the number counts of high-redshift X- 
ray quasars (§ 3.3). We subsequently compute the mass 
accreted into BHs during luminous quasar phases (§ 4) 
and compare this to the local density in BHs (§ 5), which 
we compute using the measured velocity dispersion func- 
tion of Sloan Digital Sky Survey (SDSS) galaxies (Sheth 
et al. 2003). In § 5 we also demonstrate that the most 
massive BHs were already in place at high redshift. We 
then show in § 6 and § 7 that while the Mbh — cr relation is 
redshift-independent, the ratio of BH mass to stellar mass 
within galaxies should increase with increasing redshift if 
both are regulated by a similar feedback process. Finally 
we compare our findings with some previous work in § 8 
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and summarize our main conclusions in § 9. 

2. FEEDBACK AND THE GROWTH OF GALACTIC 
BLACK-HOLES 

The growth of a supermassive BH is expected to be 
accelerated during a galaxy merger, when cold gas is 
driven to the center of the merger remnant (Mihos & 

Hcrnquist 1994; Hcrnquist & Mihos 1995). Yu & 
Trcmaine (2002) have found that the local mass density 
in BHs is consistent with the integrated luminosity den- 
sity of quasars if the mean radiative efficiency of the ac- 
creting material is erad ~ 8%; moreover, they concluded 
that about half the current BH mass density was accreted 
near the Eddington rate at redshift z > 2. It is natural 
to expect that a quasar shining at its limiting Eddington 
luminosity could generate a powerful galactic wind and 
eventually terminate the accretion process that feeds it. 

Indeed, Silk & Rees (1998) pointed out that winds from 
quasars could generate self-regulating outflows in the sur- 
rounding gas, if the energy in the outflow liberates as much 
as the binding energy of the gas in a dynamical time. 
For a quasar output power that scales as the Edding- 
ton limiting luminosity, this condition leads to a power- 
law scaling of the BH mass with halo velocity dispersion, 
Mbh oc cr^. This scaling was used by Haehnelt, Natarajan 
& Rees (1998) to demonstrate consistency between the 
density of quasars at z ^ 3 and the time derivative of the 
Press-Schcchtcr mass-function for a quasar lifetime that is 
comparable with the Salpeter time (the e-folding time of 
the BH mass), ~ 2.4 x 10'^(erad/0.06) yr. Here we extend 
an earlier model (Wyithe & Loeb 2002) to include feed- 
back limited BH growth within the context of hierarchical 
merging in a ACDM cosmology. We compare our model 
with the most recent data for th(^ space density of high 
redshift quasars at both optical and X-ray wavelengths. 

We assume that following a merger, a BH shines at a 
fraction rj of its Eddington luminosity (oc Mbh) return- 
ing fraction Fq of this energy^ to the galactic gas. The 
quasar will unbind the galactic gas if it supplies as much 
as the gravitational binding energy of the gas over the 
time it takes the gas to respond dynamically to this en- 
ergy injection, fdyn. The cold gas is usually located in a 
disk with a characteristic radius (A/\/2)''vir ^ 0.035rvir 
(Mo, Mao & White 1998), and hence a dynamical time''' of 
idyn ^ 0.035rvir/i'c- Here A ~ 0.05 is the spin parameter, 
Tvir is the virial radius, and Vc is the circular velocity of 
the galactic halo^ (Barkana & Loeb 2001), 



self-regulation condition is then. 
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where Mhaio is the halo mass, (5 Mbh is the mass accreted 
and LEdd,0 is the Eddington luminosity per unit mass 
(Af.g). Note that we assume (as appropriate at high red- 
shifts) that all the gas within a galactic halo cools on a 
time much shorter than the Hubble time. Our primary 
objective is to show that this condition is sufficient to de- 
scribe the quasar luminosity function at 2 > 2 (before 
groups and clusters of galaxies become prominent). Equa- 
tion (2) implies Afbh 3C '-'c inferred for galac;tic; halos 
in the local universe (Ferrarese 2002). The index of 5 is 
larger than the 4-4.5 inferred from the local Mbh — cr re- 
lation (Merritt & Ferrarese 2001; Tremaine et al. 2002), 
the difference having it origin in the observation that the 
Vc — c relation is shallower than linear (Ferrarese 2002). 
The index of 5 is also larger than presented in the recent 
work of Baes, Buyle, Hau & Dejonglie (2003), with the dif- 
ference having its origin in the uncertainty in the slope of 
the Mbh — cr relation. Feedback regulated growth implies 
that the Mbh — Vc relation is independent of redshift, 

M.,^1.9X10«M,(J§)(5^^)\ ,3, 

The redshift independence is consistent with the recent re- 
sults of Shields et al. (2002). Using quasars out to 2; ~ 3, 
they demonstrated that the Mbh — cr relation does not 
evolve with redshift. The resulting relation between BH 
and halo mass may be written as 

Mbh(Mhalo, Z) = eMhalo 
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Assuming 77 = 1 as suggested by observations of low 
and high-redshift quasars (Floyd 2003; Willott, McLure & 
Jarvis 2003), we can compute the normalization in equa- 
tion (4) from equation (2). The only adjustable parameter 
in our model if 77 ^ 1 is _Fq. If the BH growth was mostly 
complete by a redshift z ~ 1-2 (Yu & Tremaine 2002) so 
that the local BH masses reflect the limiting values at that 
redshift, then a value of = 10~^-^ is consistent with BH 
masses that are measured locally (Ferrarese 2002). This 
implies a value of i^q = 0.07. 

Modification of the simple scheme described above 
would be required if cooling of the gas heated by the out- 
flow is important so that additional energy is required 
to unbind the gas. The factor required could be as 
large as c/vc based on momentum conservation (Begel 
man 2003), with the resulting dependence Mbh cx: 
This relation is shallower than observed. The impor- 
tance of cooling can be evaluated by estimating the cool- 
ing time for virialized gas. The cooling rate for primor- 
dial gas of temperature Tvir lO^K and proton density 
1.2 X 10~^ [(1 -I- z)/3]'^cm^'^ (typical for a virialized halo) 

If the energy output of the quasar is communicated to the galaxy through a wind with an outflow velocity that is not much larger than the 
escape velocity of the halo, then momentum conservation leads to the same condition as equation (2). 

® We assume that gas can be driven towards the black hole during the merger. Such accretion is communicated over the characteristic dynamical 
time r/v (as appropriate for a radial orbit) rather than 2T!r/v (as appropriate for a circular orbit). 

® One inaccuracy in our modeling is that we assume dark matter halos to be singular isothermal spheres, while an NFW profile gives a lower 
Vc at the disk radius than at r-vir- We ignore this difference because galaxies appear to have nearly flat rotation curves in practice. 
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where ^(z) is close to 

Ac = IStt^ + 82rf - 39rf2, and d = - 1 (see Barkana 
& Loeb 2001 for more details). The value of Fq should be 
of order unity if the gas traps much of the quasar output 
energy without radiating it away (Begelman 2003). The 
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is A ~ 8.7 X 10-30 [(1 + z)/if (Tvir/106K)i/2erg cm-^ s'^ 
(Barkana & Loeb 2001). We have adopted a cooling rate 
that increases as \/T above T^i^ = lO^K assuming that 
it is dominated by brcmsstrahlung at higher tempera- 
tures. The thermal energy per unit volume in this gas 
is ~ 5 X 10-" [(1 + ^)/3]3 {T^.Jl{fK)evg crxT^. The 
cooling time for virialized gas is therefore 
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where the virial temperature is given by (Barkana & Loeb 
2001), 



Tvir = 2.2 X 10'' 
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3.1. The optical luminosity Junction 

We model the quasar luminosity function assuming that 
quasar activity follows major galaxy mergers (mass ratios 
greater than 3) during which BHs coalesce (Kauffmann & 
Haehnelt 2000). A linear relation between BH mass and 
halo mass (Mbh oc v^) would be a natural consequence of 
BH growth that is dominated by coalescence (Haehnelt et 
al. 1998). However, a nonlinear relation must result from 
significant mass accretion during the active quasar phase 
as found by Yu & Tremaine (2002). 

Our approach is based on the formalism described in 
Wyithe & Loeb (2002). Each quasar is assumed to shine 
at its Eddington luminosity over a lifetime tq. The quasar 
light curve may therefore be written as 



/(i) = ^|^e(i,-t). 



(8) 



The cooling time should be compared with the character- 
istic lifetime of the quasar, which our model relates to the 
dynamical time of the galactic disk. 



t,,„ = 0.035^ = 10^ [C(.)]-i/2 
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where LEdd,B = 5.73 x lO"^Mbhi0 is the Eddington lu- 
minosity in the B-band of a BH with mass Mbh in solar 
masses. The quasar luminosity function (comoving num- 
ber density of quasars per B-band luminosity Lb) is then 
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dMbh / dAMbh j 
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jiTf — "ii" , A — -rrr are the mass function and 
merger rate for BHs respectively, is the cosmic time at 
redshift z, and 5{x) is the Dirac delta function. Integrating 
over Mbh we find 
'^{Lb,z) = 
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Thus we expect the halo cooling time to be longer than 
tdyn in the range of redshifts and virial temperatures rele- 
vant for bright quasars. As soon as the gas is expelled into 
the galactic halo and is heated beyond its virial tempera- 
ture, it will not be able to cool during the quasar activity. 
Hence we assume that accretion is halted as soon as the 
quasar supplies the entire galactic gas with more than its 
binding energy. 

3. MODELING THE QUASAR LUMINOSITY FUNCTION 

The peak in the quasar luminosity function at ^; ~ 2 

coincides with the transition of the characteristic mass 
scale of collapsed halos through the mass range of inter- 
est for the Mbh — Mhaio relation. In the present-day uni- 
verse, the characteristic 1-cr mass scale is ~ 10^^-10^'' Mq. 
However galaxies reflect the characteristic mass scale of 
^ at z ^ 2, and combine to make groups 

and clusters rather than super-galaxies today. It may be 
that galaxies of this characteristic size do not form due to 
prevention of cooling through conduction or the removal 
of cold gas from halos (e.g. Benson et al. 2003), or because 
the stars in giant ellipticals have already formed and ex- 
pelled the surrounding gas at z ~ 6 (Loeb & Peebles 2003) . 
The excursion set formalism describes the rate of galaxy 
mergers near the characteristic mass at times earlier than 
z ^ 2, but describes the formation rate of groups and clus- 
ters near the characteristic mass at later times. The evo- 
lution of the bright end of the quasar luminosity function 
is therefore only tractable to analytic descriptions at early 
times (z > 2), while the faint end should be tractable at all 
redshifts. Since the observed quasar luminosity function 
implies that around half the mass in BHs was accreted 
before z ~ 2, the Mbh — Mhaio relation should be evi- 
dent in the bright end of the quasar luminosity function 
at z; > 2 — 3. We now demonstrate that this is indeed the 
case. 

^ Our model does not account for the population of obscured quasars (Fabian 2003). We implicitly assume that the powerful quasars detected 
at high redshifts are capable of removing any surrounding gaseous torus that may obscure them, since they release sufficient energy to unbind 
the gas within the entire host galaxy. 



where a = 1/(1 + z) is the scale factor, Zform is the for- 
mation redshift of the quasar, and we have used the re- 
lation between black- hole and halo mass. Evaluation of 
this expression requires the number of mergers of halos of 
mass between AM and AM -I- dAM with halos of mass 
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(Lacey & Cole 1993), as well as the 

co-moving density of halos with mass M, (Press & 

Schechter 1974). 

Assuming tq <C H^^{z), we then obtain a—Omin = j^tq, 
yielding the luminosity function'' 

/■O.SMhalo 3 
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Mhalo-AMhalo 

where we have used Mbh = L^dd ,b/(5.73 X IO^LqMqI). 
We impose a lower cutoff on the mass of halos that can 
contain a quasar. Following reionization, the cutoff mass 



4 




10 12 



10 12 
'o^ioC-b/I®) 




10 12 
logioCLe/lo) 



10 12 





1 ' 1 
z~8.0 ■ 




z~io.o ; 


A , 


1 \ i\ 1 


10 


12 14 



Fig. 1. — Comparison of the observed and model luminosity functions. The data points at z < 4 are summarized in Pei (1995), while 
the light lines show the double power-law fit to the 2dF quasar luminosity function (Boyle et al. 2000). At z ~ 4.3 and z ~ 6.0 the data 
is from Fan et al. (2001a;2001b;2003). The grey regions show the l-u range of logarithmic slope ([—2.25, —3.75] at z ~ 4.3 and [—1.6, —3.1] 
at z ~ 6), and the vertical bars show the uncertainty in the normalization. The open circles show data points converted from the X-ray 
luminosity function (Barger et al. 2003) of low luminosity quasars using the median quasar spectral energy distribution. In each panel the 
vertical dashed lines correspond to the Eddington luminosities of BHs bracketing the observed range of the Mbh-^'c relation, and the vertical 
dotted line corresponds to a BH in a 10^^'® galaxy. 



corresponds to the Jeans mass in an ionized intergalactic 
medium (IGM), i.e. a virial temperature of ~ 2.5 x lO^K 
(Barkana & Loeb 2001). Equation (2) implies that the 

quasar lifetime should be identified with the dynamical 
time of a galactic disk, tq ~ idyn, where tdyn is given in 
equation (7). 

We show a comparison between the observed and theo- 
retical luminosity functions in Fig. 1. The model success- 
fully describes the known features of the optical quasar 
luminosity function at 2 > 2. We also show on these 
plots the range of quasar luminosities corresponding to 
BH masses observed locally to define the Mbh — Afhaio re- 
lation (dashed lines). The black holes powering quasars at 
z ^ 1 are at the upper end of the mass range of the BHs 
observed locally. As discussed above, the Mbh — -Mhaio re- 
lation is obeyed locally by galaxies, but not by groups or 
clusters of galaxies inside of which the hot gas is unable 
to feed the BH. Our model luminosity function therefore 
over-predicts the number of bright quasars at low redshift. 
However, our luminosity function correctly describes the 
density of quasars at ^; ~ 1 below the break in the luminos- 
ity function. This luminosity regime corresponds to most 
of the BH mass range in the locally observed Mbh — Mhaio 
relation. In order to match the bright end of the observed 
galaxy luminosity function, the infall of gas into galaxies 
with Vc > 500 kms~^ must be suppressed (Kauffmann & 
Haehnelt 2000). This scale indeed corresponds to the cut- 
off in the Vc distribution of galaxies (Sheth et al. 2003). 
The dotted lines show the luminosity of a BH within a 
galaxy of mass Mhaio ~ 10^^'^ corresponding to a circular 



velocity of ~ 500km s~^ at z = 0. This scale is close to the 
break in the quasar luminosity function at low-redshifts. 

Near the peak of quasar activity at redshifts of 2: ~ 2-3 
our model predicts a quasar lifetime of idyn ^ 1 — 3 x 10^ 
years. This result naturally explains observational esti- 
mates of quasar lifetimes (e.g Steidel et al. 2002; Yu & 
Tremaine 2002; Haiman & Hui 2001; Martini & Wein- 
berg 2001). It is important to note that our model sug- 
gests that the observed lifetime is typical of quasars near 
the peak of the luminosity function, but that quasars at 
higher redshift have shorter lifetimes. 

At redshifts below the peak of quasar activity our model 
successfully describes the luminosity function at luminosi- 
ties fainter than the break (the agreement of the slope is 
particularly good). Boyle et al. (2000) found that the lumi- 
nosity function at 2; < 2.5 is well described by a luminosity 
evolution. However inspection of Fig. 1 suggests that the 
evolution of quasars does not follow a luminosity evolution 
at high redshifts. In particular, our model suggests that 
rather than being an intrinsic feature of the quasars them- 
selves, the break in the luminosity function is present only 
at low redshift and is a product of the inability of gas to 
cool inside massive dark matter halos. We also see that 
the faint end slope (reproduced at low redshift) evolves, 
becoming steeper at high redshift. 

We emphasize that the redshift dependence of Mbh oc 
(1 + z)^/^ in equation (4) is necessary to explain the den- 
sity of bright quasars at z ~ 6. If the characteristic feed- 
back time was constant, so that Mbh oc {\.-\- z) as is the 
case for stars (see §3), then the sharp decline in the Press- 
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Fig. 2. — Upper panels; Comparison of the observed and model X-ray luminosity functions. In each panel the vertical dashed lines corre- 
spond to the Eddington luminosities of BHs bracketing the observed range of the Mbh-t^c relation. Lower panels: evolution of the luminosity 
function with redshift in different luminosity ranges. The filled data points are from the luminosity function of Miyaji et al. (2001), while the 
open circles show results from Barger et al. (2003). 



Schechter (1974) mass function with increasing halo mass 
would have resulted in a deficit of luminous quasars that 
cannot be compensated for by arbitrarily long quasar life- 
times. For example we find that the predicted density of 
quasars at z ~ 2-3 can be made consistent with obser- 
vations by using a quasar lifetime of 10* years, but that 
the quasar abundance at 2; ~ 6 is underestimated by ~ 4 
orders of magnitude. 

The success of our model is particularly impressive, 
given the fact that its only free parameter (riFq) is di- 
rectly set by local observations of BH masses, and that 
it assumes only that BH growth in a galactic disk is self 
regulating. In the following section, we demonstrate that 
this success extends into the X-ray regime. 

3.2. The X-ray luminosity function 

Haiman & Loeb (1999) noted that the luminosity func- 
tion of quasars in the X-ray band at z ~ 3 can be re- 
lated to the B-band luminosity function at that redshift if 
the median quasar spectral energy distribution is universal 
(Elvis 1994). They then extrapolated their model to pre- 
dict the number counts of z > 5 quasars and found that 
~ 50 should be detected in a Chandra deep field. Barger 
et al. (2002) have since found one z > 5 quasar in a circle 
of radius 6' at a flux level Fx > 2 x 10" erg s'^ The 
new constraints provided by these very high redshift X-ray 
quasars (as well as optical high redshift quasars) motivate 
us to revisit the issue of modeling the X-ray luminosity 
function. 



The number counts of X-ray quasars provide strong 
additional constraints on models of the luminosity func- 
tion. At high luminosities, wide field surveys by the 
ROSAT satellite provide the X-ray luminosity function in 
the same luminosity range as the optical luminosity func- 
tion. Therefore at z 3, the comparison of the model 
optical and X-ray luminosity functions with the respective 
data sets tests the applicability of the median quasar spec- 
trum adopted for the modeling (Haiman & Loeb 1999). 
However the deep X-ray surveys now becoming available 
(Cowie et al. 2002; Barger et al. 2003) probe bolometric 
quasar luminosities at z ~ 3 — 5 that are ~ 2 orders of 
magnitude fainter than are probed in the optical. This 
large dynamic range in luminosities provides a long lever 
arm for testing models of quasar evolution. 

We follow Haiman & Loeb (1999) and construct the X- 
ray luminosity function in analogy with the B-band lu- 
minosity function using the median quasar spectrum of 
Elvis (1994). The universality of the spectrum is sup- 
ported by the latest X-ray data on high redshift quasars 
(Brandt et al. 2002; VignaH et al. 2003). In order to fa- 
cilitate comparison with observed luminosity functions we 
compute the luminosity L2-8 in an (intrinsic) X-ray band 
(e.g. 2-8keV) by integrating over the median spectrum. 
In analogy with equation (12) we find 
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Fig. 3. — Number counts of quasars in a circle of radius 6' and redshifts above 3, 5 and 7 (top to bottom) as a function of the flux limit. 
The number counts from Baxger et al. (2003) at z > 3 and z > b with flux > 2 x 10~^^ erg cm~^ are shown for comparison. 
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Here Z/2-8,0 = f2kev ^^•^^^'^>®^ where vLj.^q is the spec- 
tral energy distribution of a solar mass BH radiating at its 
Eddington limit and = hv is the photon energy. 

In the upper two panels of Fig. 2 we show the predicted 
X-ray luminosity function in the 0.5keV-2keV band (solid 
lines), and compare it with the observed luminosity func- 
tion derived by Miyaji et al. (2001) from ROSAT data 
(solid dots). Luminosity functions (solid lines) are shown 
at 2; ~ 1 and z ^ 3.5. The bolometric luminosities of 
the observed quasars are comparable with those of bright 
optical quasars. We also show results (open circles) for 
lower huninosity X-ray AGN (Cowie et al. 2002; Barger 
et al. 2003). These luminosities refer to the 2-8keV range 
and should be compared to the dotted luminosity func- 
tion. Several features of the comparison are similar to 
those of the optical luminosity function. The agreement 
is excellent at redshifts above the peak in quasar activ- 
ity. However there is again an overproduction of bright 
quasars at 2; < 2, although the luminosity function below 
the break is reproduced. The interpretation of these fea- 
tures was discussed in the previous section. As noted by 
Haiman & Loeb (1999) the success of the; modcil in produc- 
ing both the optical and X-ray quasar luminosity functions 
at 2; > 2 indicates that the median quasar spectral energy 
distribution is also a good representation at high redshifts. 
To help with the orientation between the optical (Haiman 
& Loeb 1998) and X-ray luminosity functions, we have 
placed vertical dashed lines at luminosities corresponding 
to the limits of the observed range of locally observed BH 
masses. 

In the lower two panels of Fig. 2 we show the predicted 
density of X-ray quasars (solid lines) of different luminosi- 
ties as a function of redshift in the 2keV-8keV band. These 
are compared with the observed densities derived by Cowie 
et al. (2002) and Barger et al. (2003) at low luminosities 
(open circles) and Miyaji et al. (2001) at higher luminosi- 
ties (solid dots). At 10^^-10^^ erg s~^, the model under 
predicts the quasar density at z < 1 . However at luminosi- 
ties of 10^^-10^^ erg s~^, the observed redshift evolution 
is well described by our model. 

We also compute the number counts of X-ray quasars. 
Barger et al. (2003) and Cowie et al. (2002) have recently 



presented number counts of faint, high redshift quasars 
in the Chandra deep field north. The observations were 
made in the 0.4-6keV (observed) band and the minimum 
flux level was 2 x 10"""^^ erg s~^ cm~^. We compute the co- 
moving density of quasars with observed fluxes above this 
level in this band as a function of redshift by finding the lu- 
minosity fimction in the band 0.4(1 + z) 6(1 + z)keV. This 
density is shown as the dashed line in the lower left panel 
of Fig. 2, along with the density at 2; ~ 5.7 corresponding 
to the observation of one quasar in the redshift range 5-6.5 
within a circle of radius 6' (open triangle). The model is 
in agreement with the observation of one source at 2; > 5. 

3.3. Number counts of high redshift X-ray quasars 

We can also compute the number counts of high redshift 
X-ray quasars in the Chandra deep field north. Fig. 3 
shows the number of quasars expected from the model 
with redshifts above 3, 5 and 7, and within a circle of 
radius 6' as a function of the fiux limit. These are com- 
pared with the number counts from Barger et al. (2003) 
at z > 3 and z > 5. We find that ~ 1-2 quasars with 
z > 5 and flux > 2 x 10"^^ erg s~^ cm~^ are expected 
in the Chandra deep field north compared with the one 
observed. Similarly ~ 10 quasars with z > 3 and flux 
> 2 X 10~^^ erg s~^ cm~^ are expected compared with 
9 observed. This reproduction of the high redshift, low 
luminosity number counts is an impressive success of our 
model. For example, the Haiman & Loeb (1999) model, 
which postulated a linear relation of Mbh oc Mhaioj fails 
this test. 

Motivated by this success we can use the model to pre- 
dict the number counts of X-ray quasars at even higher 
redshifts. For example the results of Fig. 3 suggest that 
0.1 quasar at z > 7 would be found in the Chandra deep 
field (or one quasar per 10 fields). Alternatively, probing 
deeper by a factor of ^ 3 in flux should allow detection of 
one quasar with z > 7 per fleld. 

3.4. Using X-ray quasars to extend the optical luminosity 
function at high redshift to lower luminosities 

The high redshift X-ray quasars discovered in the Chan- 
dra deep field have bolometric luminosities that are 2 or- 
ders of magnitude lower than the bright optically selected 
quasars. Using the median spectrum we can translate the 
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Fig. 4. — The mass density in supermassive BHs as a function of redshift for the luminosity function computed in this paper (dark line), 
and for the 2dF luminosity function (light line; Boyle et al. 2000). 



density of X-ray quasars to points on the optical luminos- 
ity function. In the 2.5 < z < 3.5 and z ~ 4.3 panels of 
Fig. 1 wc plot the density of quasars with X-ray luminosi- 
ties between lO''^ - lO*"' erg s^^ and lO''^ - lO*" erg s-\ 
and 3 < z < 5 in units of density per B-band luminos- 
ity (open circles) for comparison with the model optical 
quasar luminosity function. In the z ~ 6 panel we plot 
the density corresponding to one quasar per natural loga- 
rithm of luminosity in the redshift range 5 < z < 6.5. The 
agreement is again excellent since the comparison simply 
presents the same data in a different way. However the 
comparison demonstrates agreement between model and 
data over a large range of luminosities at high redshift. In 
particular, it extends the comparison between model and 
data to the lower logarithmic half of the observed range of 
BH masses. At z ^ 4.3 the model produces the normal- 
ization and slope of the SDSS luminosity function, then 
becomes less steep at lower luminosities, and intersects 
the point computed from the X-ray observations, some 
2.5 orders of magnitude away in luminosity (and therefore 
BH mass) and 2.5 orders of magnitude away in density 
per logarithm of luminosity. Similarly, the model satisfies 
the normalization of the luminosity function at z ^ 6 for 
both the low and high luminosity quasars as well as the 
constraint on the slope of the bright end of the luminosity 
function. The model also reproduces the suggestion from 
the data that the density of quasars changes more rapidly 
at higher luminosities between z ~ 4.3 and z ~ 6. 

4. ACCRETION OF BLACK-HOLE MASS DURING LUMINOUS 

PHASES 

The comoving mass density accumulated by quasar BHs 

by a redshift z is 

' dL^^^>{LB,z')^, (13) 

eradC-^ dz 

where Lboi is the bolometric luminosity of a quasar with 
B-band luminosity Lq , and Crad is the efficiency of conver- 
sion of accreted mass to radiation. We take e^ad = 0.06, 
appropriate for the last stable orbit of a thin disk accre- 
tion onto a Schwarzschild BH. In Fig. 4 we show the mass 
density as a function of redshift for the luminosity func- 
tion computed in this paper (dark line), and for the 2dF 
luminosity function (light line; Boyle et al. 2000) . We find 
that pbh ~ 2 X IO^MqMpc"'^ has been accreted by z ~ 1.5, 
similar to estimations of the local BH density (e.g. Yu & 



Pbh(^) = / dz' 



Tremaine 2002; Aller & Richstone 2002). Our luminosity 
function model erroneously indicates that growth contin- 
ues at z < 1.5 due to the overestimate of massive galax- 
ies at low redshift. As discussed in the analysis of Yu & 
Tremaine (2002), we find that an appreciable fraction of 
the present mass density in supermassive BHs is not ac- 
creted until after z ~ 3. Yu & Tremaine (2002) also find 
that ~ 90% of the present-day BH mass density is accreted 
before z ~ 1.5. Indeed, our luminosity function suggests 
that the accreted mass-density at z ~ 1.5 is comparable 
to the local vahie, and in agreement with the estimate 
from the 2dF luminosity function. We reiterate the point 
that by z ^ 2 tdyn is becoming comparable to the Salpcter 
time. Within the scheme outlined in this paper, that is the 
reason why most quasar growth occurs during this period. 
We therefore suggest that it is the feedback regulated ac- 
cretion around the time of the peak in quasar activity that 
shapes the Mbh — Mhaio relation. 

5. THE LOCAL MASS FUNCTION OF MASSIVE 
BLACK-HOLES 

The local density of supermassive BHs can be com- 
puted by combining the Mbh — cr relation with the velocity 
function of galaxies (e.g. Yu & Tremaine 2002; Aller & 
Richstone 2002). Until recently the velocity function had 
to be computed through combination of a galaxy lumi- 
nosity function with the Faber-Jackson (1976) and Tully- 
Fisher (1977) relations. A more reliable representation is 
now possible using the measured velocity dispersion func- 
tion of early type galaxies. Sheth et al. (2003) presented 
the measured velocity dispersion function for early type 
galaxies (which contain the most massive BHs) , as well as 
a model for the velocity dispersion function of late-type 
galaxies. They suggested an analytic fit of the form 



(t){cr)da = (/)*{ — 



exp [-{a/a^)'^] da 
r(a//3) ^ a ' 



(14) 



where (p^, is the number density of galaxies and ct^, is 
a characteristic velocity dispersion. Sheth et al. (2003) 
found that the parameters u*, a and (3 are strongly cor- 
related with one another. We take 0* = (2.0 ± 0.1) x 
10-3Mpc-3^ a = 6.5 ± 1.0, (3 = (14.75/a)°-8 and cr* = 
16ir(a//3)/r[(a + l)//3] km s'^. From Fig. 6 in Sheth et 
al. (2003) we model the velocity dispersion function in the 
regime below a 200 km s~^ (where it is dominated by 
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late-type galaxies and measurements are not provided by 
SDSS) using the functional form 



da 



<i>{a)da = logio(e) x ioi-73-i.niog,„(a) ^_ 



(15) 



To compute the mass function of BHs, we combine (/)(cr) 
with the Mbh — cf relation 



Mbh = Mbh. 



200 



( 1 V 

V2OO km s-i/ 



(16) 



.200 km s- 

Mcrritt & Ferrarese (2001) quote Mbh,200 = (1.48±0.24) x 
1O*M0 and x = 4.65 ± 0.48. We computed a large num- 
ber of realizations assuming that the parameters a, 
Mbh. 200 a-nd X a-re distributed as Gaussian within their 
quoted uncertainties. We then find the mean and variance 
of the logarithm of the resulting set of mass functions. 

In Fig. 5 we show the range (± twice the variance) in 
the computed mass function of BHs (solid area). This area 
is truncated below the region corresponding to the mea- 
sured minimum density and the maximum measured ve- 
locity bin. The grey lines show the extrapolation to lower 
densities. The derived BH mass function differs substan- 
tially from that described in Aller & Richstone (2002). In 
particular, the cutoff at the high mass end is not as sharp 
in the present case. The reason is that the observed ve- 
locity function does not agree with that derived from the 
galaxy luminosity function and the Faber- Jackson relation 
(for a discussion, see Sheth et al. 2003). 

We have computed the local mass-density in BHs by 
integrating the BH mass-function. We find a local BH 
density of pbh = (2.3^1;^) x lO^MgMpc-^, where the 
quoted error is the density computed from mass-functions 
at the boundaries of the range shown in Fig. 5. Us- 
ing the alternative values of Mbh,200 = 1-14 x io8 i3±o o6 
and X = 4.02 ± 0.32 (Tremaine et al. (2002) we find 
Pbh ^ (2.5li;^) X lO^MgMpc-^ These estimates of BH 
density agree with the estimate of mass accreted during 
luminous quasar phases as described by the 2dF luminos- 
ity function (§ 4). 

The local BH mass function may be compared with the 
redshift dependent mass function derived by combining the 
Press-Schechter (1974) mass function with the Mbh— Mhaio 
relation (4). We have plotted the resulting mass function 
at 2: = 1, 3 and 6 in Fig. 5. We find that the local density 
of BHs with masses below Mbh ^ 10*^Mq were in place by 
z ~ 1 — 3, coinciding with the peak in the quasar activity. 
More massive BHs were in place at even higher redshifts. 
In particular, we find that the most massive BHs known 
(Mbh ~ 3 X IO^Mq) may have been in place by 2; ~ 6. 
This finding is in agreement with the conclusion of Locb 
& Peebles (2003) that stars in giant ellipticals have already 
formed and expelled the surrounding gas at 2 ~ 6, but here 
the conclusion is derived from a completely independent 
argument. 

5.1. The nearest very massive black-holes 

Consider a BH of mass Mbh ^ 3 x 1O^M0. We can es- 
timate the co-moving density of black-holes directly from 
the observed quasar luminosity function and our estimate 
of quasar lifetime. At the peak of quasar activity, quasars 
powered by Afbh ^ 3 x 10*^ M0 BHs had a comoving density 
of \E'(> L) ~ 50Gpc~^. However fhubbie/^dyn ~ 2 X 10^, so 



this corresponds to the peak of BH growth, the BH den- 
sity should match the value observed today. The density 
implies that the nearest BH of ^ 3 x IO^A/q should be 

at a distance dbh (47r/3 x 10'')^^^ Gpc ~ 30Mpc which 
is comparable to the distance of M87, a galaxy known to 
possess a BH of this mass (Ford et al. 1994). 

What is the most massive BH that can be detected dy- 
namically in a local galaxy redshift survey? The Sloan Dig- 
ital Sky Survey probes a volume of ~ IGpc^ out to a dis- 
tance 30 times that of M87. At the peak of quasar activ- 
ity, the density of the brightest quasars implies that there 
should be ~ 100 BHs with masses of 3 x 10^°Mq per Gpc^, 
the nearest of which will be at a distance rfbh ~ 130Mpc, 
or ~ 5 times the distance to M87. The radius of gravita- 
tional influence of the BH scales as Mbh/i'c oc M^^^. We 
find that for the nearest 3 x IO^Mq and 3 x lO^^Mg BHs, 
the angular radius of influence should be similar. Thus the 
dynamical signature of ~ 3 x 10^°Mq BHs on their stellar 
host should be detectable. 

The luminosities of the brightest quasars at z ~ 2-3 
imply BH masses of ~ 10^*^ M©, even if they accrete at 
their Eddington rate. The masses could be lower if the 
emission is beamed, though this does not appear to be 
the case for high redshift quasars (e.g. Willott, McLure & 
Jarvis 2003). The implied density of ~ 10^°Mq BHs is in 
excess of the densities obtained by convolving the galaxy 
velocity function (Sheth et al. 2003) with the Mbh — cr re- 
lation (Merritt & Ferrarese 2001; Tremaine et al. 2002). 
This is shown by the dot-dashed lines in Fig. 5. Here we 
have estimated the density of black-holes as implied by the 
quasar luminosity function at 2; ~ 2.3 and z ^ 3 according 
to 

dMbh if-MMbh^^' ^ ^ 

where we have used dMbh = dL/5.73 x 10^. The pres- 
ence of these BHs at 2; ~ 2 — 3 may imply that these very 
massive BHs will be found in the central galaxies of large 
clusters of galaxies. 

6. REGULATION OF STAR FORMATION 

Wc may derive the limiting stellar mass of a galaxy 
in view of the same self-regulation principle (Dekel & 
Woo 2002). Suppose that stars form with an efficiency 
/* out of the gas that collapses and cools within a dark- 
matter halo, and that a fraction Fsn of each supernova 
energy output, -Esn, heats the galactic gas mechanically. 
The mechanical feedback will halt the star formation once 
the cumulative energy returned to the gas by supernovae 
equals its binding energy (assuming negligible radiative 
losses for a sufficiently rapid starburst). Hence the limit- 
ing stellar mass is set by the condition. 



M^, 1 ^6 2 

-ESN-PSN = -C/b = -7:;— Mhalot'c' 



(18) 



where wsn is the mass in stars (in Mq) per supernova ex- 
plosion. Equation (18) implies that the total mass in stars, 
M* = {fMb/ ^771) Mhs^io, scales as 

M,^A4i^ja^)]Hl + z). (19) 

At a fixed redshift, the star formation efficiency scales as 
2/3 



that the comoving density of BHs was ^ lO^Gpc-*. As A 

ocv^(x Mjj^jq. Thus smaller galaxies are less efficient at 
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Fig. 5. — The mass function of supcrmassivc BHs. The grey region shows the mass function estimated from the velocity function of Sheth 
et al. (2003) and the M^h — c relation of Merritt & Ferrarese (2001). The lower bound corresponds to the lower limit in density for the 
observed velocity function (o"^ > 10^®'^). while the grey lines show the extrapolation to lower densities. We also show the mass function 
computed at 2; = 1, 3 and 6 from the f'rcss-Schcchter (1974) halo mass function and equation (4), as well as the mass function at z ~ 2.35 
and z ~ 3 implied by the observed density of quasars and a quasar lifetime t(jyn (dot-dashed lines). 



forming stars, but a galaxy of fixed mass is more efficient 
at higlier redshift. A Scalo (1986) mass function of stars 
lias «;sN = I26M0 per supernova and -EgN = 10^^ ergs, 
and so we find that a mass in stars of = 3 x IO^^Mq 
and a velocity of Vc ^ 175km s^-'^ [the typical value ob- 
served locally (Bell & De Jong 2001)] implies i^sN ~ 0.5 
and fi, = 0.07. The star formation efficiency in the most 
massive spiral galaxies is therefore expected to be ~ 10%. 

7. REDSHIFT DEPENDENCE OF THE MAGORRIAN ET AL. 

(1998) RELATION 

In the preceding sections we have found that the mass of 
BHs scales as while the star formation efficiency scales 
as v"^. This is because the luminosity of a quasar is lim- 
ited by the Eddington value and is proportional to the BH 
mass, while the energy release rate of supernovae is pro- 
portional to the total stellar mass. Hence we find that the 
ratio of BH to stellar mass, 

^^[^{z)]Hl + z)i, (20) 

is independent of mass as observed (Magorrian et al. 

1998). However, we find that the relation evolves with 
redshift, and that the BH comprises a larger fraction of 
the stellar mass at a higher redshift. This effect has been 
seen in a sample of gravitationally lensed quasars (with 
lensed host galaxies); Rix et al. (1999) found that host 
galaxies are much fainter at z ^ 2 when compared with 
lower redshift galaxies containing quasars of comparable 
luminosity. Our results are also consistent with the semi- 
analytic predictions of Kauffmann & Haehnelt (2000). 

8. COMPARISON WITH PREVIOUS WORK 

Since the work by Haehnelt, Natarajan & Rees (1998) 
there have been several attempts to describe the quasar 
luminosity function within hierarchical merging scenar- 
ios through semi-analytic modeling. Kauffmann & 
Haehnelt (2000) assumed that the Mbh-c relation was set 
by the amount of gas that cools and becomes available 
for star formation. They found that they could explain 
the decrease in quasar activity at low redshift based on 
the decrease in merger rates, the decrease in the gas avail- 
able to fuel BHs, and the assumption that gas accrets more 



slowly at late times. For an accretion timescale that scales 
as (l-|-z)^'^/^, these considerations led to quasar hosts that 
are a factor of 10 less massive a,t z ^ 2 than in the local 
universe. Volonteri, Haardt & Madau (2002) followed the 
evolution of the BH population through a merger-tree al- 
gorithm. They modeled the luminosity function by fixing 
the Mbh — cr relation from observations, and forcing the 
quasar lifetime following a merger to be equal to the time 
for accretion of the mass to match the Mbh — relation at 
an efficiency of 10%. 

Di Matteo, Croft, Springel & Hernquist (2003) noted 
that the above studies either assume the Mbh — rela- 
tion from observations, or derive it based on additional 
assumptions. They suggested that a linear relationship 
between the gas in bulges and Mbh results in the local 
A^bh — cr relation, the drop in quasar density at low red- 
shift, and the regulation of BH growth through supernovae 
feedback. However their model over-predicts the density 
of high-redshift quasars. In their conclusion, Di Matteo et 
al. (2003) state that it seems unlikely that BH accretion 
affects gas dynamics in galaxies unless feedback from the 
quasar itself is important. On the other hand, the impor- 
tance of quasar feedback has recently been highlighted by 
Begelman (2003). 

In the previous sections, we have demonstrated that 
a scenario in which quasar feedback limits BH growth 
through regulation of the quasar lifetime naturally repro- 
duces the quasar luminosity function at high redshifts. For 
the first time we have compared the model to both the 
^; ~ 6 SDSS luminosity function and the X-ray luminosity 
function, allowing comparison over ~ 3 orders of magni- 
tude in flux and 3 orders of magnitude in quasar density 
per logarithmic luminosity interval at z > 5. The success 
of the model signals several important implications. First, 
around the peak of quasar activity at z ~ 2.5 the quasar 
lifetime, interpreted here as the dynamical time of the host 
galactic disk, is ^ 10^ years [see Eq. (7)]. This provides 
a natural explanation for observational estimates of the 
lifetime (e.g. Steidel et al. 2002; Yu & Tremaine 2002). 
Moreover, the dynamical time of the disk becomes com- 
parable to the Salpeter growth time at this redshift, so 
that most of the mass in a given BH was added during 
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this epoch. Thus we suggest that the growth of BHs dur- 
ing luminous quasar phases around this time, regulated by 
feedback from the quasar itself, set the slope and ampli- 
tude of the the Mbh — Vc relation observed today. 

9. CONCLUSIONS 

We have shown that if BHs grow in the centers of galax- 
ies until they unbind the galactic gas that feeds them over 
a dynamical time, then the locally observed relation be- 
tween BH mass and halo velocity dispersion (Fcrrarcsc 
2002) is recovered both in slope and normalization. We 
include this feedback regulation of BH growth in a model 
of the quasar luminosity function and find that the ob- 
served optical quasar luminosity function at 2; > 2 as well 
as at luminosities below the break at lower redshifts, are 
reproduced for hierarchical merging of galactic halos in 
the standard ACDM cosmology. Moreover, using the me- 
dian quasar spectral energy distribution we showed that 
the same model yields good agreement with the X-ray 
quasar luminosity function at z > 2 (and the luminosity 
function below the break at lower redshifts). The model 
suggests that ~ 10 faint quasars at ^ > 3 and ~ 1 at 
z > 5 should have been found in the Chandra deep field 
north at a flux above 2 x 10~^^ erg cm~^, in agreement 
with the observed number counts. Overall, the model de- 
scribes observations of the high redshift luminosity func- 
tion over ~ 3 orders of magnitude in bolometric lumi- 
nosity and ^ 3 orders of magnitude in comoving density 
per logarithm of luminosity. This success is impressive 
given the fact that the model is based on a simple al- 
gorithm with only one adjustable parameter, namely the 



quasar feedback power in Eddington units rjEq, and that 
this quantity is set by local observations of BH masses. 
Our model associates the inferred quasar lifetime of ^ 10^ 
years (Steidel et al. 2002) with the dynamical time of a 
galactic disk during the peak of quasar population {z ^ 2- 
3). Moreover, since the quasar lifetime is comparable to 
the Salpeter time during this epoch, the model suggests 
that the slope and normalization of the locally observed 
Mbh — Afhaio relation resulted from feedback regulated ac- 
cretion around this epoch. We have computed the local 
mass-function of BHs by combining the Mbh — cr relation 
with the measured velocity dispersion function of Sheth 
et al. (2003). Comparison of the local BH mass function 
with the mass function inferred from the feedback rela- 
tion and halo mass-function implies that the most massive 
galaxies already hosted their present-day BHs by 2: 3-6 
(see also Loeb & Peebles 2003). Finally, we showed that 
if supernova feedback regulates star formation in a similar 
way, then the Magorrian et al. (1998) relation between 
the BH and stellar masses of galactic spheroids should be 
redshift dependent [see Eq. (20)] , as observed for the hosts 
of gravitationally lensed quasars (Rix et al. 1999). 
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